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Abstract: We study how to search for photon-photon scattering in vacuum at present 
petawatt laser facilities such as HERCULES, and test Quantum Electrodynamics and non- 
standard models like Born-Infeld theory or scenarios involving minicharged particles or 
axion-like bosons. First, we compute the phase shift that is produced when an ultra-intense 
laser beam crosses a low power beam, in the case of arbitrary polarisations. This result 
is then used in order to design a complete test of all the parameters appearing in the low 
energy effective photonic Lagrangian. In fact, we propose a set of experiments that can be 
performed at HERCULES, eventually allowing either to detect photon-photon scattering as 
due to new physics, or to set new limits on the relevant parameters, improving by several 
orders of magnitude the current constraints obtained recently by PVLAS collaboration. 
We also describe a multi-cross optical mechanism that can further enhance the sensitivity, 
enabling HERCULES to detect photon-photon scattering even at a rate as small as that 
predicted by QED. Finally, we discuss how these results can be improved at future exawatt 
facilities such as ELI, thus providing a new class of precision tests of the Standard Model 
and beyond. 
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1. Introduction 

Photon-Photon Scattering (PPS) in vacuum is a still unconfirmed prediction of both Quan- 
tum Electrodynamics (QED) [1] and non-standard models such as Born Infeld theory [2, 3]. 
Additional contributions to the process can also appear in new physics scenarios involving 
minicharged [4] or axion-like [5] particles. However, all the experiments that have been 
performed by now could only be used to set upper limits on the photon-photon cross sec- 
tion <7 77 . The best constraints were obtained recently by PVLAS collaboration [6], and are 
still seven orders of magnitude above the QED prediction for a 77 . In the last few years, 
there has been an increasing interest in studying the possibility to detect PPS at future 
facilities, using two possible strategies. On one hand, the cross section for the process will 
be maximum at a possible future photon-photon collider [7], based on a free electron laser 
producing two beams of photons in the MeV range. A second approach will be to perform 
experiments at optical wavelengths, and compensate the smaller cross sections with a very 
high density and/or a long path of interaction of the colliding photons [8, 9, 10, 11]. Most 
of these proposals require the construction of new facilities, that will eventually be available 
in the future, such as a free electron laser and/or an exawatt laser. Two exceptions are 
Refs. [10, 11], that discussed the possibility of performing experiments at present facilities 
to improve the PVLAS limit on the photon-photon cross section. However, even in these 
cases the predicted sensitivity was found to be sufficient to detect PPS of QED origin only 
at future facilities, such as ELI [12] or VIRGO+ [13], respectively 
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Here, we will propose a set of experiments that can already be performed at present 
petawatt laser facilities, such as HERCULES [14, 15]. First, we perform a new theoretical 
computation to get a quantitative expression for the phase shift that is produced by PPS 
when two orthogonally polarised beams cross each other. This result turns out to provide 
a test for PPS of QED origin which is significantly more sensitive than our previous pro- 
posal [9, 10], in which the two crossing beams had the same polarisation. Moreover, it can 
be used in combination with our previous result to provide a full test of all the parameters 
appearing in the low energy effective Lagrangian describing the photons in non-standard 
models, such as Born Infold theory and scenarios involving minicharged (MCP) or axion- 
like (ALP) particles. In fact, taking into account the precision that can be achieved in 
the measurement of optical phase shifts or ellipticities [16, 17, 6], we propose a set of ex- 
periments that will allow either to detect PPS at HERCULES, or to set new limits on 
the relevant parameters, improving by several orders of magnitude the current constraints 
obtained by the PVLAS collaboration. We then propose a multi-cross optical mechanism 
that can further improve the sensitivity of this set of experiments, eventually enabling 
HERCULES to detect PPS as predicted by QED. Finally, we discuss how these results can 
be improved at future exawatt facilities such as ELI, thus providing a new class of precision 
tests of the Standard Model and beyond. 



2. The effective Lagrangian for the electromagnetic fields in QED and in 
non-standard models 

We will consider the case of photon energies well below the threshold for the production of 
electron-positron pairs, and assume an effective Lagrangian for the electromagnetic fields 
E and B of the form 

C = C Q + i L Cl + -^ T G\ (2.1) 

being Co = y ^E 2 — c 2 B 2 ^ the Lagrangian density of the linear theory, Q = eoc(E • 
B) and eo and c the dielectric constant and the speed of light in vacuum, respectively. 
The additional, non-linear terms, that appear multiplying the parameters £l and £r in 
equation (2.1), are the only two Lorentz-covariant terms that can be formed with the 
electromagnetic fields at the lowest order above Cq. Therefore, they will appear as the first 
correction to the linear evolution both in QED and in non-standard models. 

In fact, in QED photons can interact with each other through the interchange of virtual 
charged particles running in loop box diagrams [1]. Besides other interesting effects [18], 
such an interaction leads to the Euler-Heisenberg effective Lagrangian density [19], that 
coincides with equation (2.1) with the identification 5,^° = S^ ED = £, being 



8a 2 h 3 



3 



« = ii^- 6 - 7xl0 "T- (2 - 2) 

On the other hand, in Born- Infold theory [2], we would obtain the relation 
4£|^/7 [3], in general without a definite prediction for the numerical value. 
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The presence of a minicharged (or milli-charged) particle (MCP) [4] would provide an 
additional contribution analogous to that from the electron-positron box diagram. If the 
new MCP are spin 1/2 fermions, and assuming that their mass m e is still larger than the 
energy of the photons (the eV scale in optical experiments), we would obtain 

A# CT = A# CT =(i^) 4 f, (2.3) 

where e is the ratio of the charge of the particle with respect to the electron charge. The 
existing laboratory bounds in this regime is e < 8 x 10~ 5 [4]. Taking masses above the 
eV scale, in order to apply the effective Lagrangian approach, this limits can be read as 
^^MCP < O (io 6 ,f). As we shall see in the next sections, this constraint has been improved 
by PVLAS collaboration, and can be further strengthened by the experiments that we 
propose in the present paper. Of course, in this case there are already stronger limits, 
e < 1CT 15 , from astrophysical and cosmological observations [4]. A larger contribution 
might be obtained if the MCP are lighter than the energy scale of the photons (the eV 
scale in the present paper). However, this case would deserve a different treatment which 
goes beyond the purposes of the present work, since it cannot be described simply by an 
effective Lagrangian of the form of equation (2.1). 

Similar considerations apply if the new MCP is a spinless boson. Assuming again 
that its mass m e is still larger than the energy of the photons (the eV scale in optical 
experiments), and using the results of Ref. [20], we would obtain 

16 \ m e J 

and 

= (2 - 5) 

On the other hand, if the MCP is a spin 1 boson, the contribution to the effective La- 
grangian would be larger, as computed using the result of Refs. [20]. We obtain 

mcp 1 = 261(-iM 4 (26) 
16 \ m e J 

and 



a^ icpi = — fi^Ve (2.7) 



28 V m e . 

Let us now discuss the case of an axion-like particle [5]. This can be a Light Pseu- 
doscalar Boson or a Light Scalar Boson, depending on the coupling with the photons, 
that is described in the Lagrangian density by the terms Cp = — \fhcgp§pQ and Cs = 
—VKcgs^s^-G-, respectively. We can find the leading contribution to the effective La- 
grangian when the photon energy is much smaller than the m<j> scale, that can be cast in 
the form of equation (2.1) with an additional contribution given by 

*T = ^4 (2.8) 
7cm J 
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and A£l = 0, in the case of pseudoscalars, or 

Ml = p\ (2.9) 

and A£r = 0, in the case of scalars. On the other hand, for m$ > leV, the Cristal Ball [5] 
laboratory limit g P < 4.2 x 10 _3 GeV _1 gives the contraint gp/(m$c 2 ) <4x 10 6 GeV~ 2 , 
which can be converted in the limit A£r < 2.2 x 10~ 25 m 3 /J. This constraint has been 
improved recently by the PVLAS consideration [6], as we shall see in the next section. 
Again, the astrophysical limits gp < 2.7 x 10 -9 GeV _1 , valid for m$ < IKeV [5], is still 
much more stringent than any laboratory bound. 

Similar considerations apply for scalar boson, for which the best laboratory constraints 
are also those that were recently set by PVLAS, that we will review in the next section. 
We also recall that our approximations do not apply for masses smaller than the order 
of the energy of the colliding photons. In this case, the computation of is more 
complicated, and the production of real axions has also to be taken into account. The 
latter can be expected to produce dichroism, just as in the presence of a constant external 
magnetic field [21, 6], and combined with the measurement of ellipticity may allow for 
a determination of both m$ and gs,p- However, this case lies beyond the scope of the 
present paper, that uses a phenomenological approach that can be applied to any theory 
that goes beyond the Standard Model, in the energy regime in which it only implies a 
different contribution to equation (2.1), as parametized by artitrary £l and £t- Expressing 
the electromagnetic fields in terms of the four-component gauge field = (A , A) as 
B = V A A and E = — cVA° — ^A, this gives the equations of motion as the variational 
derivatives ST/SA^ = 0, where r = / Cd A x is the effective action. Such equations are 
similar to the modified nonlinear Maxwell's equations that have been obtained in Ref. [22], 
the only difference being the distinction between £^ and 



3. Present constraints 

The current limits on PPS in vacuum have been obtained recently by the PVLAS collab- 
oration [6] by searching evidence of birefringence of the vacuum in a uniform magnetic 
field background [21]. Their negative result was used to set the current constraints on the 
parameters appearing in equation (2.1). With our notation, their 95 % C.L. limit reads 

|7 ^- 4al < 3.2 x 1(T 26 ^. (3.1) 

Assuming £l = £r = £ exp as in QED, their results can be translated in the limit ^ exp < 
3.2 x lCT 26 m 3 /J, which is 4.6 x 10 3 times higher than the QED value of equation (2.2) (7 
orders of magnitude for the cross section <r 77 ). 

Note however that PVLAS experiment was only sensitive to the combination |7£t — 4£j,| 
of the parameters. In particular, this quantity vanishes when £p = therefore PVLAS 
experiment is unable to set any constraint on a pure Born Infeld theory 
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4. Approximate solution for the scattering of orthogonally polarised beams 

In Ref. [9], we have studied the scattering of two counter-propagating waves that are 
polarised in the same direction, and we have found that the effect of PPS was to produce 
a phase shift in each wave, which was proportional to £l multiplied by the intensity of the 
other wave. That result was obtained by an analytical, variational approximation, and was 
shown to agree with a numerical solution of the full non-linear equations, that was also 
obtained in the second of Refs. [9]. 

Here, we will apply a similar variational method to find a solution for the problem of the 
scattering of two orthogonally polarised counter-propagating waves, one of which represents 
an ultra-high power beam. Let the low power and the high power waves be polarised in 
the x and y directions respectively, so that their linear (free) evolution (neglecting photon- 
photon scattering) would be A x (t, z) hn = ao cos(kz-tot) and A y (t, z) hn = Acos(kz+ut+ip) , 
where u = kc and we allow for an initial phase difference (p. Their energy density, when 
each of the two waves is taken alone, would be p x = e w 2 a 2 /2 and p y = p = € lu 2 A 2 /2, 
respectively. Hereafter, we will assume that p y » p x . When these two waves are made 
to scatter, they will affect each other due to PPS, as described by the non-linear terms in 
equation (2.1). First, we note that the assumption of no dependence on x and y of the 
fields guarantees that the condition At = A z = is maintained by the non-linear evolution, 
since we have checked that in this case the equations 5T/5A t = and ST/5A Z = are 
automatically satisfied, independently on the values of A x (t, z) and A y (t, z). Therefore, in 
the absence of x and y dependence, the components A t and A z with will not be generated 
if they are not present from the beginning. Second, we note that the non-linear effect is 
driven by the very small parameters and This justifies a perturbatively-motivated 
variational approach, similar to that introduced in Refs. [9]. We then need to chose a good 
ansatz for the fields A x (t, z) and A y (t, z). In principle, each of the two components can get 
a transmitted wave contribution, propagating along the same direction as the original wave, 
and a reflected wave propagating in the opposite direction. In a perturbative approach, we 
can compute these different effects separately and then sum them up. Therefore, we will 
first neglect the reflected waves, and use the following ansatz: 

A x = a{z) cos(kz — cut) + (3(z) sin(kz — cot), (4.1) 
A y = A cos(kz + ut + if). 

Here, we have neglected the effect of the low power wave on the high power wave, taking 
into account that such an effect is expected to be proportional to the energy density of the 
low power beam. This expectation, inspired by our previous work [9], will be confirmed 
by the result that we will obtain below. 

We now substitute the ansatz (4.1) in the Lagrangian (2.1), and average out the fast 
variation in z over distances of the order 2ir/k, assuming that the envelop functions a(z) 
and (3{z) will show a much slower variation, as we will verify a posteriori. We then compute 
the variational equations 5T/5a = and 5T /5(3 = 0, keeping the lowest order terms in the 
expansion parameter £t and neglecting the higher order space derivatives (d/dz) n of a(z) 
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and /3(z), as compared with k n . After a long but straightforward algebra, we find the 
following equations: 

P , (z)+x T a(z)=0, (4.2) 
a'(z) - xtP(z) = 0, 

where xt = 7eQC 2 A 2 k 3 ^T/^- Assuming the initial condition a(0) = ao, (3(0) = 0, in such a 
way that the corrected solution coincides initially with that of the linear problem, we find 
then a(z) = aocos(xrz), (3(z) = — ao sin(xr-?)- After substituting in equation (4.1), we 
finally get the following variational solution: 

A x = aocos(k z + xt z — uit), (4.3) 
A y = A cos(kz + u)t + ip) . 

In other words, taking into account that eoA 2 uj 2 /2 ~ p is the energy density of the 
high power wave, we find that after a crossing distance Az, the phase of the orthogonal, 
low power wave is shifted by an amount A<pT = XtAz = 7^TpkAz. Note that this result 
is independent of the initial phase difference (p between the two crossing waves. 

Let us now introduce the possibility that a reflected wave is generated in the component 
A x , as described by the ansatz 

A x = ao cos(kz + xtz — cot) + j(z) cos(kz + tot) + S(z) sin(fcz + uit), (4.4) 
A y = A cos(kz + cot + ip) 

After repeating the same procedure as above, we get the following variational solution 

7 (z) = ^ [cos(147t P £t) - 1] , (4.5) 
S( Z ) = ^ sin (i4vr^ T ). 

Now, the quantity lAnp^x can be estimated for the petawatt laser HERCULES [15] that 
we will consider below for our proposals of experiments. In this case, the peak intensity 
is / ~ 2 x \0 22 WcmT 2 [15], corresponding to an energy density p ~ 6.7 x 10 17 JmT 2 . 
Taking into account the PVLAS limit £ exp < 3 x 10~ 26 m 3 /J, and assuming that it can be 
applied to £t at least roughly (see also figure 2 in the last section), we find that £tP ^ 
2 x 1CP 8 for the product giving the importance of the nonlinear QED effects. As a result, 
the variational solution for 7 and 5 implies that j(z) ~ for all the practical purposes, 
and S(z) ~ 7aop^Tkz/2. Taking k ~ 7.8 x 10 6 m -1 as for the wavelengths of e.g. the 
HERCULES laser (A = 800ram), we obtain that the \5(z)/a \ < O.hz/lm. Now even in the 
multi-cross configuration that will be discussed below the crossing length will be smaller 
than the centimetre scale, so that 6(z) will be smaller than the low power amplitude Qo at 
least by two orders of magnitude. For this reason, it will be neglected. 

Finally, let us introduce the possible reflected wave in the component A y , as described 
by the ansatz 

A x = Qo cos(kz + xtz — Lot), (4.6) 
A y = Acos(kz + cot + ip) + rj(z) cos(kz — ut) + a(z) sin(kz — tut). 
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By repeating the same kind of computations and arguments as above, we find the following 
solution 



A x = ao cos(kz + xtz — tot), 
A y = A cos(kz + ojt + ip) + i]o cos(kz + \lz — cot) , 



(4.7) 



where xl = 2eoc 2 ^l 2 /c 3 ^L = 4^p/c, 77(0) = 770 and we assume that <r(0) = 0. We then see 
that the counter-propagating wave in the y polarisation only exists if it is present from the 
beginning, and that it gets a phase shift A^l = xlAz which is equal to that obtained in 
Ref. [9], as could be expected. 

equation (4.7) implies that, after crossing a counter-propagating, linearly polarised 
ultra-intense laser pulse, an ordinary laser pulse is phase shifted both in the polarisations 
parallel and orthogonal to that of the high power beam. The corresponding phase shifts 
are 



where I = pc \s the intensity of the high power beam and r = Az/c is its time length. 
If we assume £l = £r as in QED, we see that A4>t is more sensitive by a factor 7/4 
than A(pL to the effect of PPS. This is already an improvement with respect to Ref. [9]. 
Moreover, the dependence of equations (4.8) on both parameters £l and £77 will permit 
a full analysis of the effective Lagrangian (2.1), distinguishing between QED and other 
models such as Born Infeld theories. Finally, we note that (4.8) also implies that the 
high power pulse behaves like a birefringent medium, producing a relative phase shift 
A(j)b = AipT — A4>l = (7£r — 4£x)/A;t between the transverse an parallel polarisations of 
the low power beam. 

5. Proposal of experiments 

We will now discuss how the result of equations (4.8) can be used to search PPS in vacuum 
by measuring phase shifts and ellipticities. In fact, Ref. [16, 17] provides a technique 
that allows for the measurement of phase shifts as small as 10~ 8 rad, which is the noise 
limit [17]. This precision, that holds for ultra-short laser pulses [16, 17], applies then to our 
A<pL and A<pT- A similar sensitivity can be obtained for the measurement of the ellipticity 
induced by A<^ corresponding to birefringence. In particular, in the same experiment that 
we have cited above [6], the PVLAS collaboration was able to resolve the corresponding 
A(j)b with a statistical error 05 = 1.1 x 10~ 8 rad, thus allowing them to set a 95% C.L. 
experimental limit Acfii, < 2.8 x 10~ 8 rad. Hereafter, to be definite and for simplicity, we 
will use this same numerical value, 2.8 x 10~ 8 rad, for the sensitivity in the measurement 
of A(f>L and AcpT, taking into account that the actual experimental precision will be close 
to this choice [16, 17]. 

We can now propose a set of three experiments: 



A^l = ^ L pkAz = 4Z L IkT, 
A<pT = T^rpkAz = 7£, T IkT, 



(4.8) 
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1) A linearly polarised low power laser pulse is divided by a beam splitter in two 
branches. One of them propagates freely in vacuum, while the other crosses a contra- 
propagating ultra-high power laser pulse polarised in the same direction. The phase shift 
suffered by the low power pulse as a consequence of PPS is then measured by comparing 
with the pulse that has propagated freely, using the technique described in Ref. [16, 17]. 
Due to equation (4.8), this configuration can be used to measure the parameter £x = 4 ^ ^ , 
where we have introduced a gain factor F that corresponds to a multi-cross configuration 
as discussed below. This experiment will then allow either to detect PPS by measuring a 
non- vanishing £l, or to set an upper limit on this parameter as 

2.8 x 1(T 8 , . 

< ~wmr- <"> 

2) The configuration is the same as in case 1), except that now the high power beam is 
polarised in a direction orthogonal to that of the low power pulse. Due to equation (4.8), 
this setup can be used to obtain the parameter £r = -/f tkr ^Y measuring the phase shift 
A(()t- This will allow either to detect PPS or to set the upper limit 

2.8 x 10" 8 /r n . 

3) The low power beam polarisation has two components, one parallel and another 
orthogonal to that of the contra-propagating high power pulse. Ellipticity measurements 
can then be used to deduce the difference of the phase shifts A<f>T — A^l = Acftb, allowing 
to determine the combination 7 ^ T ~ 4 ^ = , thus allowing either to detect PPS, or to 
set the upper limit 

|7fr-4fr| 2.8 x KT 8 

3 < 3FIkr ■ [b - 6) 
The combination of these three experiments will permit a complete exploration of the 
parameter space. Actually, it is easy to see that if £l an d £r have the same sign, as in 
QED and Born Infeld theories, experiment 3) is less sensitive than the combination of the 
others and can be discarded without losing significant information. On the other hand, if 
£l and have opposite sign, experiment 3) is the most sensitive one, although even in 
this case the other two measurements would be useful for a full determination of both £l 
and £t- 

From equations (5.1), (5.2) and (5.3), we see that the sensitivity depends on the com- 
bination Ikr of the experimental parameters of the ultra-high power laser beam, and on 
a gain factor F that will be discussed later. Therefore, the most favourable experimental 
configuration will be that allowing for the maximum value of the product Ikr. As far as 
we know, the highest value achieved at present facilities is that of HERCULES laser [15], 
that reaches peak intensities I = 2 x 10 22 W/ cm 2 , for a time length r = 3 x 10 _14 s and 
wavelength A = 2ir/k = 8.1 x 10~ 7 m. This gives Ikr = 4.7 x 10 19 J/m 3 . Even in the 
absence of any gain factor (F = 1), such a facility will be able to resolve £l and £t as small 
as 1.5 x 10 _28 m 3 / J and 8.6 x 10 _29 m 3 /J, respectively, thus allowing either to detect PPS 
of non QED origin, or set limits on the parameters that are more than two orders of mag- 
nitude (5 order of magnitude in the cross section) more stringent than the current PVLAS 
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limits, in addition to the fact that they constrain the full parameter space, including the 
case of Born-Infeld theories, that were unconstrained by PVLAS. 

A significant improvement will be obtained in the near future at ELI [12], that in its 
first stage will achieve peak intensities I ~ W 25 W/cm 2 , for a time length r ~ 10~ 14 s and 
wavelength A = 2-ir/k ~ 8 x 10~ 7 m. This gives Ikr ~ 8 x 10 21 J/m 3 . Even in the absence 
of any gain factor (F = 1), such a facility will be able to resolve and £t as small as 
9 x 10~ 31 m 3 /J and 5 x 10 _31 m 3 /J, respectively. In particular, ELI would allow for the 
detection of PPS of QED origin and for measuring the parameter £ with two figures. 

6. Improving the sensitivity with multiple crossing 

The sensitivity of our proposed experiments can be enhanced by making the two beams 
cross each other several times, using a kind of wave guide consisting of two parallel series 
of parabolic mirrors as shown in figure 1. An advantage of using parabolic mirrors is that, 
in the paraxial approximation, they do not generate aberrations in the beam. We assume 
that the laser pulses are localised at a distance R half a way to the path leading to the 
next mirror. To be concrete, we will also assume that at the crossing points the high 
power laser is focused to the diameter d ~ 0.8/um and intensity / ~ 2 x \0 22 Wcm~ 2 of the 
HERCULES beam. The time duration r = 30/s implies that the pulse length cr ~ 9/wn 
along the direction of propagation is approximately an order of magnitude greater than 
its transversal width, therefore the two beams must cross forming an angle 9 close to it 
in order to maximize their superposiposition. These requirements may be achieved using 
plasma mirrors [23], that can work e.g. at the intensity /mirror — 2 x 10 19 W/cm 2 with a 
reflection coefficient r ~ 0.98 [23]. In this case, the high power beam at the mirror should 
have a diameter equal to d^J I// m i rr0 r — 25/xm. In order to avoid diffractive distortions, we 
will use parabolic mirrors of, say, a double diameter, d m i rr0 r — 50/xm. On the other hand, 
the two planes are assumed to be at a distance ~ 2R, with R much larger than d m i rr0 r, say 
R = 5cm, in such a way that 9 — tt = 2 aiccos(d m i rroi / R) — tt ~ — 2 x 10~ 3 , so that 9 is very 
close to 7r. A serious technological challenge to be faced will be the very precise alignment 
of the mirrors, since any uncertainty in the direction will be multiplied by the number of 
times the beams are reflected. In principle, the orientation of the mirrors can currently be 
fixed with a precision as small as A9 ~ 10~ 8 rad [24]. After TV reflections, this will produce 
an uncertainty ~ NRA9 on the position of the spot at the focus. This uncertainty must 
be smaller than the diameter of the beams at the focus, so that TV < d/(RA9) ~ 10 3 . 

In order to compute the gain factor F of this configuration, we note that after each 
reflection the intensity is reduced by a factor r ~ 0.98. Moreover, the phase shift is 
due to the counter-propagating components of the photon momenta, p z = hks'm(9/2). 
Taking into account that k appears to the third power in the expression of the phase 
shift (or equivalently that it appears to the sixth power in the cross section [1]), the gain 
factor is then F = s'm 4 (9/2)J2n=o rn > wri ere we have included a further factor sin(0/2) 
taking into account that Az becomes cr sin(0/2) in this configuration. This result can 
also be obtained in more elegant and rigorous way by making the computation in the 
reference system in which the total momentum is zero and the two colliding photons are 
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Laser 1 



Plasma mirrors 



antiparallel. In fact, by indicating with a prime the quantities in such a system, and 
being z and y the vertical and horizontal directions in thelaboratory system of figure 
1, we have: t' = j(t — (3y/c), y' = 7(7/ — (3ct), z' = z, to' = 7(0; — (3ck y ), and k! y = 
j(k y — (3uj/c) =0 and k' z = k z , where (3 = ck y juj and 7 = \J\ — (3 2 - It is then easy to see 
that the phase k' z z' — uj'if + 2eQj?(k' z f'^L/S.z l ', when translated to the laboratory system, 
gives k z z + k y y — ojt + A</>, with A(f> = sin 4 (6 , /2)A0e =7r . A similar result can be obtained 
in the case of orthogonally polarized waves. 

Taking N = 1000 in the expression of F 
that we have obtained above, we can find a 
limiting value -F max ~ 50, that in principle 
can be achieved with present technology. As 
a result, the measurement of the phase shifts 
in the experiments 1) and 2) that we proposed 
above with such a gain factor will be able to 
resolve £j, and £t as small as 3.0 x 10~ 30 m 3 /J 
and 1.7 x 10~ 30 m 3 /J respectively, thus allow- 
ing to detect PPS as predicted by QED, or 
find a signal of non-standard physics. Note 
that the combination of the two experiments 
will also be able to test Born Infeld theory 
and scenarios involving MCPs or axion-like 
particles, taking into account the discussion 
of section 2. 




Laser 2 



Plasma mirrors 



Figure 1: Proposed setup for multiple crossing 
of the two scattering laser pulses. 



7. Conclusions 

We have computed the phase shifts affecting a low power laser beam that crosses an high 
power laser pulse with general transverse polarisation, and proposed a set of experiments 
to completely determine the parameter space of the effective Lagrangian that describes 
PPS well below the threshold for the creation of electron-positron pairs. 

Our results are summarised in figure 2, showing the 95% C.L. exclusion regions that 
can be obtained with this set of experiments at the present facility HERCULES, without 
or with multi-crossing, as compared to the current constraint by PVLAS. The predictions 
of QED and Born-Infeld theory are explicitly indicated. Additional contributions from 
minicharged particles, or axion-like scalar or pseudoscalar bosons, can sum with them and 
produce a different point in the £l and £r plane, as discussed in section 2. 

Even with the single-crossing version (F = 1), figure 2 shows how the PVLAS limits 
can be substantially improved at HERCULES, possibly allowing for the detection of PPS of 
non-standard origin. On the other hand, by using a multi-cross mechanism, HERCULES 
would already be able to detect PPS of QED origin. Note that the result of figure 2 
is obtained using the conservative value F = 30, corresponding to just N = 44 aligned 
mirrors, which is a more realistic assumption than the maximum value that we have found 
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above. Finally, in figure 2, we also see how the sensitivity will be improved at ELI in the 
future, thus allowing for a more precise determination of the parameters. 

We think that this proposal 
can eventually contribute to a new 
class of precision tests of QED and 
non-standard models, such as Born- 
Infeld Theory or scenarios involv- 
ing minicharged particles or axion- 
like, scalar or pseudoscalar bosons. 



8. Acknowledgements 

We thank Eduard Masso for ex- 
plaining us the contraints on axion- 
like particles, and Frank Wise and 
Guido Zavattini for clarifying the 
sensitivities that can be achieved 
in the measurements of phase shifts 
and ellipticities, respectively. We 
are also grateful to Heinrich Hora 
and Luis Plaja for helpful remarks 
and to Bernardo Adeva for use- 
ful discussions. D.T. thanks Pe- 
dro Fernandez de Cordoba and the 
whole InterTech and IMPA Groups 
for their very nice hospitality at 
the Universidad Politecnica de Va- 
lencia. This work was partially sup- 
ported by Xunta de Galicia, by contract No. PGIDIT06PXIB239155PR, and by the Gov- 
ernment of Spain, by contracts N. TIN2006-12890, FIS2007-29090-E and FIS2008-0100. 

References 

[1] V. Costantini, B. De Tollis, G. Pistoni, Nonlinear Effects in Quantum Electrodynamics, Nuovo 
Cim. 2A (1971) 733. 

[2] M. Born, On the Quantum Theory of the Electromagnetic Field, Proc. R. Soc. London 143 
(1934) 410; 

M. Born and L. Infeld, Foundations of the New Field Theory, Proc. R. Soc. London \AA 
(1934) 425. 

[3] V. I. Dcnisov , I. V. Krivchenkov, and N. V. Kravtsov, Experiment for measuring the 

post-Maxwellian parameters of nonlinear electrodynamics of vacuum with laser-interferometer 
techniques, Phys. Rev. D 69 (2004) 066008. 

[4] B. Holdom, 2 U(l)s and epsilon- charge shifts, Phys. Lett. B 166 (1986) 196; 



- 11 - 



10- 25 i 




10 -31 10-30 10 -29 10 -28 10 -27 10 -26 10 -25 



Figure 2: Exclusion plot for the search of PPS. The pa- 
rameters £l and £t are measured in units of m 3 /J. The 
diagonal line corresponds to Born Infeld theory, while the 
point is the QED prediction. The darkest region is excluded 
by the current PVLAS constraint. The next two inner re- 
gions correspond to the parts of the parameter space that 
can be probed at HERCULES with single crossing (F = 1) 
or multiple crossing (choosing F = 30), respectively. Fi- 
nally, the last inner region represents the range that can be 
reached at ELI with single crossing. 



B. Holdom, Searching for epsilon charges and newU(l), Phys. Lett. B 178 (1986) 65; 
E. Golowich, R.W. Robinett, Limits on millicharged matter from beam dump experiments, 
Phys. Rev. D 35 (1987) 391; 

M. I. Dobroliubov and A. Y. Ignatev, Millicharged particles, Phys. Rev. Lett. 65 (1990) 679; 
S. Davidson, B. Campbell and D. Bailey, Limits on particles of small electric charge, Phys. 
Rev. D 43 (1991) 541; 

S. Davidson, S. Hannestad, G. Raffclt, Updated bounds on milli-charged particles, J. High 
Energy Phys. 05 (2000) 003; 

H. Gies, J. Jaeckel, A. Ringwald, Polarized light propagating in a magnetic field as a probe for 
millicharged fermions, Phys. Rev. Lett. 97 (2006) 140402; 

M. Ahlers et al, Particle interpretation of the PVLAS data: Neutral versus charged particles, 
Phys. Rev. D 75 (2007) 035011. 

[5] E. Masso and R. Toldra, Light spinless particle coupled to photons, Phys. Rev. D 52 (1995) 
1755; 

J. A. Grifols, E. Masso and R. Toldra, Gamma rays from SN 1981 A due to pseudoscalar 
conversion, Phys. Rev. Lett. 77 (1996) 2372; 

J. W. Brockway, E. D. Carlson and G. C. Raffelt, SN 1981 A gamma-ray limits on the 
conversion of pseudoscalar s, Phys. Lett. B 383 (1996) 439; 

E. Masso and R. Toldra, New constraints on a light spinless particle coupled to photons, Phys. 
Rev. D 55 (1997) 7967; 

E. Masso, Axion and Axion-like particles, Nucl. Phys. 114 (Proc. Suppl.) (2003) 67; 

P. Arias, J. Gamboa, H. Falomir, F. Mendez, Effects of heavy bosonic excitations on QED 
vacuum, Mod. Phys. Lett. A 24 (2009) 1289. 

[6] M.Bregant et al., PVLAS, Limits on Low Energy Photon-Photon Scattering from an 
Experiment on Magnetic Vacuum Birefringence, Phys. Rev. D 78 (2008) 032006. 

[7] I. Ginzburg, G. Kotkin, V. Serbo, V. Tclnov, Production of high-energy colliding 77 and je 
beams with a high luminosity at VLEPP accelerators, Pizma ZhETF 34 (1981) 514; JETP 
Lett. 34 (1982) 491. 

[8] S. L. Adler, Photon splitting and photon dispersion in a strong magnetic field, Ann. Phys. 
(NY) 67 (1971) 599; 

Y. J. Ding and A. E. Kaplan, Nonlinear magneto- optics of vacuum: Second-harmonic 
generation, Phys. Rev. Lett. 63 (1989) 2725; 

F. Moulin and D. Bernard, Four-wave interaction in gas and vacuum: definition of a 
third- order nonlinear effective susceptibility in vacuum: chi( (3)) (vacuum), Opt. Comm. 164 
(1999) 137; 

D. Bernard et al, Search for stimulated photon-photon scattering in vacuum, Eur. Phys. J. D 
10 (2000) 141; 

E. Lundstrom et al., Using high-power lasers for detection of elastic photon-photon scattering, 
Phys. Rev. Lett. 96 (2006) 083602; 

G. Brodin, M. Marklund, and L. Stenflo, Proposal for detection of QED vacuum nonlinearities 
in Maxwell's equations by the use of waveguides Phys. Rev. Lett. 87 (2001) 171801; 

E. B. Aleksandrov, A. A. Anselm, and A. N. Moskalev, Double-refraction of vacuum in the 
field of intense laser-radiation, Sou. Phys. JETP 62 (1985) 680; 

T. Hcinzl et al., On the observation of vacuum birefringence, Opt. Comm. 267 (2006) 318; 
A. Di Piazza et al., Light diffraction by a strong standing electromagnetic wave, Phys. Rev. 
Lett. 97 (2006) 083603; 



- 12 - 



M. Marklund and P. K. Shukla, Nonlinear collective effects in photon-photon and 
photon-plasma interactions, Rev. Mod. Phys. 78 (2006) 591; 

F. Ehlotzky et at, Fundamental processes of quantum electrodynamics in laser fields of 
relativistic power, Rep. Prog. Phys. 72 (2009) 046401; 

H. Gies, Strong laser fields as a probe for fundamental physics, Eur. Phys. J. D 55 (2009) 311; 
M. Marklund and J. Lundin, Quantum vacuum experiments using high intensity lasers, Eur. 
Phys. J. D 55 (2009) 319; 

S. V. Bulanov et al, Relativistic laser-matter interaction and relativistic laboratory 
astrophysics, arXiw.0812.1^21. 

[9] A. Ferrando, H. Michinel, M. Seco, and D. Tommasini, Nonlinear phase shift from 
photon-photon scattering in vacuum, Phys. Rev. Lett. 99 (2007) 150404; 
D. Tommasini, A. Ferrando, H. Michinel, M. Seco, Detecting photon-photon scattering in 
vacuum at exawatt lasers, Phys. Rev. A 77 (2008) 042101. 

[10] D. Tommasini, A. Ferrando, H. Michinel, M. Seco, Nonlinear Propagation of Crossing 
Electromagnetic Waves in Vacuum due to Photon-Photon Scattering, FRONTIERS IN 
MODERN PLASMA PHYSICS, AIP Conference Proceedings 1061 (2008) 292. 

[11] G. Zavattini, E. Calloni, Probing for new physics and detecting non-linear vacuum QED 
effects using gravitational wave interferometer antennas, Eur. Phys. J. C 62 (2009) 459. 

[12] http://www.extreme-light-infraestructure.eu. 

[13] G. Losurdo et al., Advanced Virgo Preliminary Design, VIR089A08. 

[14] G. A. Mourou, T. Tajima, and S. V. Bulanov, Optics in the relativistic regime, Rev. Mod. 
Phys. 78 (2006) 310. 

[15] HERCULES, V. Yanovsky et al., Ultra-high intensity-300-TW laser at 0.1 Hz repetition rate, 
Optics Express 16 (2008) 2109. 

[16] I. Kang, T. Krauss, and F. Wise, Sensitive measurement of nonlinear refraction and 
two-photon absorption by spectrally resolved two-beam coupling, Opt. Lett. 22 (1997) 1077. 

[17] F. Wise, private communication. 

[18] D. Tommasini, Reality, measurement and locality in Quantum Field Theory, J. High Energy 
Phys. 07 (2002) 039; 

D. Tommasini, Photons uncertainty solves Einstein- Podolsky- Rosen paradox, Opt. Spectrosc. 
94 (2003) 741. 

[19] O. Halpern, Scattering processes produced by electrons in negative energy states, Phys. Rev. 
44 (1933) 855; 

H. Euler, The dispersion of light to light according to Dirac's theory, Ann. Physik 26 (1936) 
398; 

W. Heisenberg and H. Euler, Folgerungen aus der Diracschen Theorie des Positrons, Z. Physik 
98 (1936) 714. 

[20] S. I. Kruglov, Pair production and vacuum polarization of vector particles with electric dipole 
moments and anomalous magnetic moments, Eur. Phys. J. C 22 (2001) 89. 

[21] R.Baier and P.Breitenlohner, The vacuum refraction index in the presence of external fields, 
Nuovo Cim. B47 (1967) 117; 

S.L.Adler, Photon splitting and photon dispersion in a strong magnetic field, Ann. Phys. (NY) 
67 (1971) 559; 



- 13 - 



Z. Bialynicka-Birula and I. Bialynicka-Birula, Nonlinear effects in quantum electrodynamics. 
Photon propagation and photon splitting in an external field, Phys. Rev. D 2 (1970) 2341; 
E.Iacopini and E.Zavattini, Experimental-method to detect the vacuum birefringence induced 
by a magnetic-field, Phys. Lett. B 85 (1979) 151. 

[22] J. McKenna and P. M. Platzman, Nonlinear interaction of light in a vacuum, Phys. Rev. 129 
(1963) 2354. 

[23] H. Hora, Plasma mirror for High Constrast Picosecond laser Pulses for Fast Ignition Fusion, 
J. Phys: Conf. Series, 112 (2008) 022028. 

[24] Ke-Xun Sun, P. Lu, R. L. Bycr, Improved Grating Angular Sensor for LISA and MGRS, J. 
Phys: Conf. Series, 154 (2009) 012033. 



-14- 



